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DIRECTIONAL VARIATIONS OF THE NON-GAUSSIANITY PARAMETER /, 
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ABSTRACT 

We investigate local variations of the primordial non-Gaussianity parameter /nl in the WMAP 
data, looking for possible influence of foreground contamination in the full-sky estimate of Jnl- We 
first improve the needlet bispectrum estimate in ( Rudjord et al. 2009) on the full-sky to /jvl = 73±31 
using the KQ75 mask on the co-added V+W channel. We find no particular values of /jvl estimates 
close to the galactic plane and conclude that foregrounds are unlikely to affect the estimate of /jv l in 
the V and W bands even for the smaller KQ85 mask. In the Q band however, we find unexpectedly 
high values of /jvl in local estimates close to the galactic mask, as well as significant discrepancies 
between Q band estimates and V/W band estimates. We therefore conclude that the Q band is too 
contaminated to be used for non-Gaussianity studies even with the larger KQ75 mask. We further 
noted that the local Jnl estimates on the V+W channel are positive on all equatorial bands from the 
north to the south pole. The probability for this to happen in a universe with /jvl = is less than 
one percent. 

Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 



1. INTRODUCTION 

Although the Cosmic Microwave Background (CMB) 
fluctuations are usually assumed to follow a Gaussian 
distribution, the consensus is that this is only an ap- 
proximation valid to a certain level of precision. The 
deviations from Gaussianity seem to be small and are 
therefore difficult to estimate precisely. However, accu- 
rate knowledge of both the size and nature of such a de- 
viation would be of high value in cosmology. The search 
for such anomalies in the CMB has therefore attracted 
much attention recently. 

Most models for inflation predict the CMB fluctua- 
tion to be slightly non-Gaussian. This level of non- 
Gaussianity is meas ured by the parameter /jvl (see e.g. 
(jBartolo et al.1120041 )). Alternative inflationary scenarios 
predict different values of Jnl, thus making an accu- 
rate estimate of this parameter crucial for understand- 
ing the physics of the inflationary era. Several /n l pa- 
rameters have been considered in the literature, corre- 
sponding to different ansatz for the shape of primordial 
non-Gaussianities. In this paper we focus on local non- 
Gaussianity, parametrized by /j^£ ai (simply f^L in the 
follow ing) . 

In (Ru diord et al~l [2009L we used needlets to esti- 
mate a value of Jnl = 84 ± 40 (local type) differ- 
ent from zero at the 2a level. This was in agree- 
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ment with previous estimates (lYadav fc Wandeltl [20081: 
i Komatsu et al.l 120091: iSmith et al. I I2009t iCurto et al. I 
120081 120091: iPietrobon et al. 1 120081 12001 where values 
of /jvl deviating from zero with about 2a were also re- 
ported. Also in agreement with previous estimates, we 
found that the value of /jvl increases with smaller sky 
cuts. This could be a result of smaller error bars due to 
the increased amount of data, or it could be an indication 
of foreground residuals near the galactic plane influenc- 
ing the value of Jnl- This is our main motivation for 
estimating /jvl on smaller parts of the sky. By having 
local estimates of /jvl we can check whether the high 
value of /jvl comes from areas close to the galactic plane 
thus indicating the influence of foreground residuals, or 
whether consistent /jvl values are found in different parts 
of the sky. The second motivation for studying the di- 
rectional dependence of Jnl, is due to several reports 
of the CMB deviating from statistical isotropy in various 
wavs (|Hansen et al. ll2008UGroeneboom fc Eriksen|| 200 "' 
i Hoftuft et all 120091: lEriksen et al.l l2004t lHansen et al 
120041 iVielva et al.l 12004 iTegmark et"a!Tl2003l) . Here we 
check whether a similar asymmetry is seen in the value 

Of f N L. 

In (IRudiord et al. I I2009T) . the bis pectrum of spheri- 
cal needlets (jLan fc Marmuccill2008[ ) was used to obtain 
/jvl- The localization and uncorrelation properties of 
needlets make them a conve nient tool for stud ying local- 
ized regions on the sky (see (jBaldi et al. | |2006)). We will 
here use a similar analysis to estimate the level of non- 
Gaussianity on selected regions on the sky and in this 
way study the spati al variations of f?yr,. We therefore 
refer the reader to (jRudiord et al. I (|2009D ) for a more 
detailed description of the procedure, as well as for a 
more extensive list of references. 

This paper is organized as follows. In section [2] we 
describe the data used for the analysis, in section [3] we 
describe the non-Gaussian maps and in section|4]we out- 
line the method used for estimating /jvl- We apply the 
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procedure to WMAP data and present the results in sec- 
tion [5] before we summarize and conclude in section [6l 

2. DATA 

This analysis was performed using the foreground re- 
duced co-added V+W frequency bands of the WMAP 
five year data at HealpbQ resolution N S id e = 512. Also 
the individual Q, V and W bands were used for consis- 
tency tests. The Gaussian simulations were generated 
using the best fit power spectrum from the WMAP five 
year release. We also used the beam and noise properties 
supplied by the WMAP team. For masking out galac- 
tic foregrounds we used the KQ75 as well as the smaller 
KQ85 mask. 

3. NON-GAUSSIAN MAPS 

The non-Gaussian maps used for the analysis have 
been generated usi n g the algorit h m de veloped in 
(ILiguori et al I 120031: ILiguori et al. I I2007T ) (see also 
lElsner fc Wandelt I (|200S^ ) for recent developments). We 
briefly review the general structure of the algorithm here, 
while addressing the reader to the above mentioned pa- 
pers for further details. The CMB multipoles ai m are re- 
lated to the primordial gravitational potential $ through 
the well known formula: 




¥ $(k)Yi m (fc)A,(fc) , (1) 



where Ag(k) is the radiation transfer function and 
the potential is written in Fourier space. However, for 
the case of primordial non-Gaussianity, the primordial 
potential takes a very simple expression in real space, 
where: 

$(x) - <Mx) + /nl [$|(x) - <*i(x))] . (2) 

In the previous expression <£>l is a Gaussian ran- 
dom field, characterized by a primordial power spectrum 
P{k) = Ak n ~ 4 ; in the following we will refer to $l(x) 
as the Gaussian part of the primordial potential. The 
remaining non-Gaussian part of the potential is simply 
the square of the Gaussian part in each point (mod- 
ulo a constant term, necessary to enforce the condition 
($(x)) = 0; however it is clear that this term affects only 
the CMB monopole) 

It is then convenient to work directly in real space and 
recast formula (fT]) in the following form 

a lm = J d 3 x<S>(x)Y em (r)A e (r) , (3) 

where Ag(r) = J dkk 2 je(kr)Ae(k) are the real space 
transfer functions, je(kr) is a spherical Bessel function, 
and r is a lookback conformal distance. This formula 
suggests to structure the algorithm in the following steps 

1. Generate the Gaussian part $i of the potential in 
a box whose side is the present cosmic horizon. 

2. Square the Gaussian part point by point to get the 
non-Gaussian part. 

7 http://hcalpix.jpl.nasa.gov 



3. Expanding in spherical harmonics the Gaussian 
and non-Gaussian parts of the potential for differ- 
ent values of the radial coordinate r in the simula- 
tion box. 

4. Convolve the spherical harmonic expansions of $l 
and $nl with the radiation transfer function Ag (r) 
in order to obtain the Gaussian and non-Gaussian 
part of the multipoles of the final NG CMB sim- 
ulation. For a given choice of the non-Gaussian 
parameter /nl a CMB map is then obtained sim- 
ply through the linear combination ag m = af m + 

fwhafm (the superscripts L and NL always indi- 
cating Gaussian and non-Gaussian respectively) . 

The most difficult and time consuming part in this pro- 
cess is actually the generation of the Gaussian part of the 
potential $. The difficulty arise from the fact that we are 
working in a box of the size of the present cosmic horizon, 
that in conformal time is about 15 Gpc, but at the same 
time a cell in this box must have a size no bigger than 
20 Mpc in order to resolve the last scattering surface, 
where most of the CMB sig n al is generated . It t urns 
out (see ILiguori et all pOOl : ILiguori et al~l ([20071 ) for 
details and more explanations) that a convenient way to 
achieve this is to work directly in spherical coordinates, 
use a non uniform discretization of the simulation box 
(since no points are needed in a large region of the box 
where photons are just free streaming) and generate the 
multipoles of the expansion of $l (x) through the follow- 
ing two step approach: 

1. Generate uncorrelated radial multipoles ri£ m (r), 
gaussianly distributed and characterized by the fol- 
lowing spectrum: 

(ne imi (n)n*e 2m2 (r 2 )) = ^F^fc i ( 4 ) 

where S D is the Dirac delta function. 

2. Filter the multipoles rii m with suitable functions 
in order to produce a Gaussian random field with 
the properties of the multipole expansion of the 
primordial Gaussian potential $l . It can be shown 
that the expression of the filter functions is: 



We(r,n) = - j dkk 2 y/P*(k)jt(kr)jt(kn) , (5) 

where P$ is the primordial power spectrum, and 
the filtering operation takes the form 



$L( r ) = / dr x r\ n im (n)W e (r, n) . (6) 

In the last expression (r) are the desired quan- 
tities i.e. the multipoles of the expansion of the 
Gaussian part of the primordial potential for a 
given r. 
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4. METHOD 

The bispectrum has shown to be the most powerful tool 
for estimating Jnl- The bispectrum is zero for a Gaus- 
sian held, and any significant deviations from a zero bis- 
pectrum is therefore a non-Gaussian signal. Needlets are 
a n ew type of spherical wav elets which were introduced 
by (jNarcowich et al. 1 120061 ). For o ur analysis we use 
the b ispectrum of needlet coefficients (jLan fc M arinucci 
2008) to estimate /nl- The localization properties of 
needlets make it possible to obtain the bispectrum in 
several different regions of the sky with very little addi- 
tional costs in terms of CPU time compared to one full 
sky analysis. Although the needlet bispectrum does not 
yield optimal error bars on f NL (sec ( Smi th et al. 112009ft 
for the optimal method), the advantage is the possibil- 
ity of a fast and easy calculation of local estimates for 
consistency checks. 

The needlet coefficients are denoted f3^ k , j is frequency, 
k the direction on the sky (we will take k as the pixel 
number in the Healpix grid) and the parameter B char- 
acterizes the localization in frequency domain. Indeed, 
needlets allow for a tight control of localization in har- 
monic space and uncorrelation in pixel space; these prop- 
erties are valuable for statistical inference and are not 
shared by other w av elet constructions, plea se refer to 
(jBaldi et al. 1120061) h (IMarinucci et al. II2008T) for details 
and further references (see also (jPietrobon et al~l 120061) 
and (|Fav et al. II2008D ). For instance, for a given value 
of B a needlet coefficient only contains information on 
multipoles in the range i = [B 3 , B 3+ ]. Thus, the pa- 
rameter B controls localization in harmonic space: small 
values of B correspond to small ranges of frequencies j, 
while the reverse is true for larger B. 

The needlet bispectrum with base B may be expressed 



then simply 



pixclsGdO 
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(7) 



where aft, are the standard deviation of the needlet co- 

j k 

efficients j3? k for a Gaussian map and dfi is the region of 
the sky for which the bispectrum is calculated. 

For the estimation of / j yx we used a similar proce- 
dure as in IRudiord et al. I (|2009l ) , but with a few small 
changes. We used a higher value for B, resulting in fewer 
needlet scales. Intuitively one would expect this to give 
a poor result, since each needlet scale would cover a large 
interval in I -space. However, a higher value for B also 
gives better localization properties in pixel-space, thus 
minimizing the influence of the mask and thereby reduc- 
ing the error bars. Additionally, this greatly reduces the 
computational cost of the analysis. 

We will calculate the needlet bispectra using the 
needlet coefficients calculated on the pixels k of a iV S idc = 
512 map. We will use all the iV s ide = 512 pixels k which 
are inside an A^ide = 2 pixel as our smallest region dfi. 
We thus obtain 48 bispectra Ij 1 j 2 j 3 (p) where p denotes a 
pixel of the -/V s idc = 2 pixelization. Having obtained these 
48 bispectra it is straightforward to construct the bispec- 
trum on sky patches with different shapes and sizes. We 
see from eq. [7] that the bispectrum of a larger patch is 



JjuWa (region) = ^^^(p), 



(8) 



where the sum over p goes over the -/V S id c — 2 pixels 
within the desired region. 

In order to estimate f^L we perform a x 2 analysis for 
every map and region to be investigated 



X 2 (.fNL) = d T (/ ArL )C- 1 d(/ 7VL ) 



where the data vector is 



d = I? h U-{I hhh (f NL )) 
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(9) 



(10) 



(Ijijzfa) i s here the average first order non-Gaussian 
bispectrum obtain ed usin g non-Gaussian sim ulations 
(jLiguori et al. Il2007| ). See (jRudiord et al. 1120091 ) for de- 
tails. The corresponding covariance matrix C is evalu- 
ated by means of Monte-Carlo simulations. 

Differentiating to find the value for Jnl which gives 
the lowest x 2 yields the "Generalized Least Squares" es- 
timate: 

T 



/jVL — 



Ijl3233 ) C ^313233 



Ijl3233 ) C 1 \ ^313'- • ; 



(11) 



We estimate loc al Jnl values accordin g to the proce- 
dure described in (Rud iord et al. Il2009i ). but with one 
important difference: We calculate and save the bispec- 
tra calculated on each individual N S id e — 2 pixel for each 
simulation. Then, when the bispectrum and thereby /jvl 
is estimated for a larger region, the bispectra for the dif- 
ferent N s ide — 2 pixels are added up according to eq. [5] 
and the correlation matrix is constructed from this final 
bispectrum for each larger region. 

5. RESULTS 

First we estimated f^L on the full sky using B = 1.781 
and 11 needlet scales covering multipoles up to £ ma x 
1 500. This was the best trade-off (lowest error bars) 
between needlet coefficients for low values of B being 
more affected by the mask but having more frequencies 
j and for high values of B being less affected by the mask 
but having fewe r scales. Otherw i se we followed the same 
procedure as in IRudiord et al. 1 (|2009f ). The results are 
presented in table [TJ 



frcq. channel 


KQ75 


KQ85 


V + W 


73 ±31 


78 ±29 


V 


58 ±35 


55 ±33 


w 


74 ±37 


72 ±34 


Q 


-8 ±39 


-9 ±37 



TABLE 1 

/jVL ESTIMATES AND lcr ERROR BARS. THESE ESTIMATES ARE 
FOUND USING THE FULL CMB SKY. 



As we see from the results, the combined V + W chan- 
nel gives a 2a deviation from Gaussianity. The individual 
V and W channels are consistent (within 2a), but the Q 
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channel deviates significantly from the others, suggesting 
possible contamination by foregrounds. 

We repeat the above analysis on localized regions on 
the sphere. First the bispectra were found for an ensem- 
ble of Gaussian simulations on each of the 48 (N S id e — 2) 
Healpix patches on the sky. We then combined these 
patches in three different ways according to eq. [5J 

1. Hemispheres The larger regions were defined to 
be hemispheres. For each of the 48 directions de- 
fined by the {N S id e = 2) pixel centers, /jvl was 
estimated on a hemisphere centered on this direc- 
tion. 

2. 45° discs Same as for the hemispheres, estimat- 
ing /jvl on discs with 45° radius instead of hemi- 
spheres. 

3. Equatorial rings These regions were defined to be 
the 7 constant latitude rings in which the N a a e = 2 
healpix pixels are ordered. Each of these rings 
covers a relatively small part of the sky, and the 
error-bars are therefore large. However they are 
not overlapping, and the /jvl estimates on the dif- 
ferent rings can therefore be considered nearly in- 
dependent (except for the largest scales). 

The Jnl estimates from the 45° discs with the smaller 
KQ85 mask are shown in figure [TJ where each pixel rep- 
resents the /jvl estimate on a disc centered on the pixel. 
The lower plot shows the same map, but /jvl for each re- 
gion has been normalized by its standard deviation. We 
see that even for the smaller mask, there is no evidence 
for particularly high values of /jvl in the galactic region, 
but note that most of the values are positive. In fig- 
ure [5] we show the corresponding map of /jvl estimates 
over hemispheres while using the KQ75 mask. As an ad- 
ditional test of consistency we have produced the same 
estimates based on bispectra obtained from iV s idc = 4 
pixels instead of iV s i<ie = 2 pixels. The result is shown in 
the same figure. We can see that the two maps show the 
same structures on the sky. 

As seen from figure[T]the distribution of /jvl on the sky 
shows a dipole. For the /nl/c map of the 45° discs of 
the V+W channel (KQ85), the dipole has a maximum at 
9 = 129°, 4> = 96° with an amplitude of T/j^ = 61. We 
investigated whether such a dipolar distribution of /jvl 
was common in Gaussian simulations and found the value 
of the dipole amplitude to be in good agreement with 
simulated maps. The dipole is therefore to be expected. 
Investigation of the hemisphere results as well as for the 
KQ75 mask using both the combined V+ W channel and 
the individual Q, V and W frequency channels yielded 
similar results. 

In order to check whether the 48 estimates of /jv l were 
internally consistent, we combined the estimates to form 
a data vector d = [f^ L , f% Ll fff L ] for a \ 2 test. Sim- 
ulations were then used to find a covariance matrix con- 
taining the covariances between the 48 estimates. Then a 
X 2 test was performed on Gaussian simulations as well as 
on the WMAP data. The % 2 values of the data were com- 
pared with those of the simulations. The results, both 
for the hemisphere estimates as well as for the 45° disk 
estimates, was that the x 2 of the WMAP data was fully 




-QMS ^^^h 2.5 



Fig. 1. — The upper figure shows /jvl estimates on 45° discs 
centered on the given pixels, while the lower figure shows the same 
estimates divided by their standard deviation. The estimates were 
made on the V+W channel using the KQ85 galactic mask. 




Fig. 2. — Both the figures shows /jvl estimates on hemispheres, 
but for the upper one the estimates are performed while using 
resolution N a id e = 2, while for the lower one we have used N a i^ B = 
4. We can see that the two maps show the same structures. These 
estimates were made on the V+W channel using the KQ75 galactic 
mask. 

consistent with Gaussian simulations (within 2a). The 
local /jvl estimates are therefore internally consistent. 

As a further test we obtained estimates of /jvl on equa- 
torial rings. The motivation behind this approach was to 
uncover possible foreground contamination outside the 
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KQ75 mask. The /jvl estimates for the different rings 
are presented in table [2] We see that the estimates of 
/jvl seem higher around equator, but the errors-bars are 
also larger (due to the galactic mask in the equatorial 
region). Since none of the rings show particularly high 
values (compared to the error bars) we do not have ev- 
idence to claim that foreground residuals have an influ- 
ence of the estimates of /jvl • The internal consistency of 
the various ring estimates have been tested as described 
above and shown to be in agreement with simulations 
(within 2a). 

We also estimated /jvl on the equatorial rings of the V 
+ W map using the smaller KQ85 mask. These results 
are presented in table [2j As expected the error-bars are 
somewhat smaller then the results with the KQ75 mask, 
especially around the equatorial region. The estimates 
are consistent with the ones from the KQ75 mask. 

We then followed through with similar investigations 
of the individual Q, V and W frequency channels using 
the KQ75 mask (table [2]). The V and W channels are 
consistent (within 2a), while only the Q channel shows 
a 3cr deviation in the ring around equator exactly where 
foreground residuals would be expected. We thus sus- 
pect a possible influence of foreground contamination in 
this band. This was further checked by testing the con- 
sistency of the estimated /jvl between the bands using a 
X 2 approach. Whereas the (V-W) and (V-Q) differences 
were found to be consistent with simulations (within 2a), 
the (W-Q) difference was found to be larger than in 99% 
of the simulations. 

Also, a similar test was performed on the differences 
f NL (V) - f NL (VW) and f NL (W) - f NL (VW). We see 
in table Q] and on some rings in table [2] that the VW 
estimates seem driven by the W estimate. This is not 
seen in simulated maps and it was found that the small 
difference f NL (W) - Jnl{VW) found for the WMAP 
data is found in only 5% of the simualted maps whereas 
the difference f n l{V) - f n l{VW) for WMAP was found 
to be consistent (well within 2a) with simulations. 





V+W 


V+W 


Q 


V 


W 


ring 


KQ75 


KQ85 


KQ75 


KQ75 


KQ75 


1 


91 ±95 


93 ±95 


47 ± 118 


66 ± 106 


94 ± 111 


2 


11 ±68 


6 ±68 


-18 ±83 


1±76 


28 ±79 


3 


80 ±80 


43 ±71 


-149 ±100 


19 ±90 


13 ±93 


4 


283 ± 183 


122 ± 113 


700 ± 226 


462 ± 205 


253 ± 213 


5 


117 ±82 


128 ± 70 


-61 ± 103 


37 ±92 


122 ± 96 


6 


39 ±66 


53 ±66 


-81 ±81 


35 ± 74 


15 ±78 


7 


158 ± 93 


156 ± 93 


174 ± 114 


138 ± 104 


201 ± 108 



TABLE 2 

The /jvl estimates and la error-bars for equatorial rings 
OF THE individual Q, V, AND W frequency channels as well 

AS THE CO- ADDED V+W FOR THE MASKS KQ75 AND KQ85. 

We see that all of the rings for the V+ W channel (as 
well as the individual V and W channels) give a estimate 
Inl > 0. These rings are nearly independent (except 
for the largest scales) and for a sky with /jvl = one 
would expect that each of these have a 50% probability of 
being positive. For 7 rings one would then expect only a 
57 ~ 0.78% probability that all rings give /jvl > (this 
probability is confirmed by simulations). 



Since the lowest value estimated for the rings of the V 
+ W channel (with the KQ75 mask) is /jvl = ff, we in- 
vestigated the probability of this occurring in Gaussian 
simulations. We found that of 4000 Gaussian simula- 
tions, only 0.35% have /jvl > H in ah the rings. 

6. CONCLUSIONS 

In this paper we have used the bispectrum of needlets 
to obtain local estimates of /jvl on the WMAP five year 
data. We performed the analysis on the combined V+ W 
channel, as well as the individual Q, V and W channels, 
using multipoles up to i — 1500 and the KQ75 galactic 
cut. For the combined V+W channel we also applied the 
KQ85 mask. 

We first made a full sky analysis, resulting in a best fit 
value of /jvl = 73 ± 31 for the combined V+W channel 
using the KQ75 mask. The individual V and W chan- 
nels give consistent (within 2a) results, but the /jvl esti- 
mate of the Q channel deviates significantly, suggesting 
contamination of foregrounds. 

The estimates were then made on selected regions of 
the sky and showed how the needlet bispectrum approach 
is powerful for finding estimates of /jvl hi many different 
regions, roughly at the cost of one single full sky estimate. 
We divided the sky into smaller regions according to four 
different patterns: hemispheres, 45° disks and equatorial 
rings. In each of these schemes /jvl was estimated in 
every region. The results were compared to simulations 
using a % 2 test, and all local /jvl estimates were found 
to be internally consistent (within 2a) for the V and W 
channels. 

The local estimates of /jvl showed a dipolar distribu- 
tion of /jvl on the sphere, with a maximum at 9 = 129°, 
<p = 96°. For comp arison, the hemispher ical power asym- 
metry reported in Han sen et al. I (|2008t ) was found with 
a maximum in 6 = 107°, <fi = 226°, and is therefore not 
expected to be connected to the findings in this paper. 

Also, such a dipolar distribution was found to be com- 
mon in simulated Gaussian maps. For the equatorial 
rings we found a positive value for /jvl in every ring. 
We compared the lower estimate of the V+W channel 
{Inl > 11 for all rings) with simulations and found that 
only 0.35% of Gaussian simulations have this feature. 

For the rings we find no significant evidence of fore- 
ground contamination outside the galactic KQ75 mask. 
The results also seem to be fairly consistent (within 2a) 
between the individual V and W channels whereas the 
Q band show signs of possible foreground contamination 
in the equatorial band where /jvl is larger than zero at 
the 3tr level. This is confirmed by the fact that the dif- 
ferences in local /jvl values between the Q and W bands 
are larger than in 99% of the simulations. 

We conclude that our study shows no significant 
anisotropy in the estimates of /jvl in the CMB sky. No 
abnormal values for /jvl are found close to the equator 
except for the Q band where we suspect foregrounds to 
influence the estimate of /jvl- 
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